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The T=0 phase diagram of the planar t-J model with small doping is investigated 
by exact diagonalisations of square clusters with up to 32 sites. At half-filling, the 
O | single hole quasi-particule weight vanishes in the strong correlation limit J — ► 

O while the quasi-particle mass diverges (like 1/J). However, the spectral function 

shows weight distributed over a large energy range (~ 7t) and a pronounced struc- 
ture in momentum space which supports the composite particle picture of the hole. 
For J/t > J/t\s,2 the effective attractive force between holes leads to bound-pairs 
formation. The hole-hole and hole pair-hole pair binding energies calculated on the 
^26 x V26 cluster indicates that in the range J/t £ [0.16,0.45] a phase of sepa- 
rate pairs of d-wave internal symmetry is stable. On the other hand, the d-wave 
£^ pair spectral function exhibits a quasi-particle peak down to very small J/t ratios 

suggesting that J/t\s,2 could in fact be as small as 0.05. Above J/t\s,4 ~ 0.45 the 
2-holes pairs bind into J r holes pairs. An abrupt change of the orbital symmetry 
ON of the 4-hole droplet is also observed above J/t ~ 2.7. In the parameter range, 

J/t\B,2 < J/t < J IAb,a, Bose condensation of the individual pairs is expected to 
^3 lead to d-wave superconductivity. 

S PACS numbers: 7/,.72.-h, 71.27. +a, 71.55.-i 

1 

5 1. INTRODUCTION 

o 

^ 1.1. Strong Coupling Hamiltonian 



Soon after the discovery of superconductivity in the cuprate materials P.W. 
Anderson 1 and others suggested that this phenomena was intimately linked to mag- 
netism. The t-J model, 

H = ^xycJ^Cy^ + ^JxySx-Sy, (1) 

x,y,<7 x,y 

defined on a two-dimensional (2D) square lattice, emerged as one of the simplest 
model to describe the low energy physics. The hole hopping matrix elements t X y 
and the magnetic exchange couplings J X y (of magnitude t=l and J respectively) 
connect only nearest neighbor sites. Hamiltonian (1) is a strong-coupling version of 
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the well-known Hubbard model. It can also be derived from the more complicated 
CuC>2 Hamiltonian in some relevant parameter regime. 2 In the following we shall 
present the state-of-the-art of the numerical work on the 2D t-J model. 

Table 1 



Sizes of the full and reduced Hilbert spaces (corresponding to the ground-state) for 
various N-sites t-J clusters with Nh holes. 



N h 1 \[~N x VW 


Full Hilbert space 


Symmetry group 


Reduced Hilbert space 


1/4x4 
1/a/26 x ^26 
1/a/32 x ^32 
4/^26 x ^26 


102 960 
135 207 800 

9 617 286 240 

10 546 208 400 


Tie 

T 32 <X> C 4v 

T 26 ®C A ®I 2 


6 435 
5 200 300 
37 596 701 
50 717 244 



1.2. Exact Diagonalisation Technique 

The Exact Diagonalisation (ED) technique is based on the fact that lattice 
Hamiltonians like (1) defined on N-sites clusters (with N a spin-u electrons) are finite 
sparse matrices. Obvious conservation laws allow one to work at fixed magnetization 
#1 — N[ (chosen to be minimum) and number of holes Nh = A — Aj — Aj. Although 



Fig. 1. Single hole and d-wave pair QP weights vs J (from Refs. 4, 10 and 12). 
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the full Hilbert space of size Af = N\/ (N-^lN^lNhl) grows exponentially fast with 
the number of sites (see Table 1), the number of matrix elements increases only 
as N x Af instead of Af 2 . The Lanczos algorithm which requires a small number 
of recursive applications of the Hamiltonian matrix can then take advantage of its 
sparseness. To minimize finite size effects the clusters are chosen with the symmetry 
of the lattice (square s/N x s/N toruses with periodic boundary conditions). As 
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shown in Table 1 use of the symmetry groups (translation group Tjv, point groups 
C4 or C$ v , and spin inversion 1 2) allow one to conveniently split the matrix into 
several blocks to be handled separately. ED have been extensively applied to the 
study of both static and dynamical correlations. 3 

2. SINGLE HOLE IN THE ANTIFERROMAGNET 

At half-filling the system develops long range antiferromagnetic (AF) order 
with low energy spin-wave excitations. Upon doping it becomes a metal with a 
Drude weight proportional to the hole doping. 3 Many studies were devoted to the 

Fig. 2. Single hole spectral function at momenta (0,0), (71", 7r) and (tt/2, tt/2) on 32- 
and 26-sites clusters (from Ref. 10). 
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single hole problem 3 which captures most of the physics of the low doping limit. In 
this case, it becomes possible to perform finite size scaling. 4,5 These studies show 
consistently (i) a quasi-particle (QP) band at low energy (see Fig. 1) and (ii) a 
substantial fraction of the spectral weight of the single hole propagator at higher 
energy. The minimum of the QP dispersion occurs at the point S (q = (7r/2,7r/2)) 
and is almost degenerate with the X point (q = (71", 0)) 4 ~ 7 in surprisingly good 
agreement with the simple Hartree-Fock treatment. The bandwidth is reduced to 
~ 2.2 J. 4,5 In other words the QP mass is enhanced by the correlations. 

Contrary to one-dimension (ID), a string-like force exists in 2D between the 
hole charge and the hole spin centers of mass. 8 However, finite size scaling 4 shows 
that the QP weight desappears as ~ J 6 when J — ► (see Fig. 1) suggesting that 
the hole eventually breaks apart in this limit. 9 On can then interpret the QP as 
a bound state between two different objects, a light one (with a mass ~ and 
a heavy one (with a mass ~ 1/=/) in analogy with the ID case. 9 The confining 
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string-like force vanishes in the J — ► limit where the two components become 
eventually free. This picture implies then naturally that simultaneously (i) the QP 
mass diverges when J — ► 3 ' 4 and (ii) the spectral function (see Fig. 2) exhibits 
weight over an energy scale of several t as well as some well-defined structure in 
k-space. 5, 10 Moreover, this is also consistent with the fact that the optical mass 
(extracted from the Drude peak) oc 1/t and the QP mass oc 1/J are essentially 
different. 11,5 

Fig. 3. A4 for various symmetries of the 4-hole GS. The short-dashed and long- 
dashed lines correspond to the exact large-J extrapolation and to the estimation 
based on the exact Heisenberg magnetic bond energy respectively (from Ref. 18). 
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3. BOUND PAIRS FORMATION 

For increasing J, the QP pole becomes stronger. A simple broken-bond counting 
argument shows that, at large J/t, two holes injected in the AF can minimize the 
local magnetic energy by forming a bound state. Finite size scaling have shown 12 
that this picture is actually correct even down to J = 3b;i — 0-2 (see also Fig. 4). 
The hole-hole pair has an d x 2_ y 2 orbital symmetry although it exhibits dominant 
hole-hole correlations at an intermediate distance of \/2 ie when the holes stay across 
the diagonal of a plaquette on the same sublattice. 13 The two-hole bound-state also 
exhibits flux quantization in a ring in units of hc/'Ie. 11 

When pairs are formed in this parameter regime they are likely to Bose conden- 
sate into a superconducting state unless some sort of clustering occurs. 14 However, 
early ED 15 as well as high temperature expansions 16 seem to consistently find that 
phase separation (PS) occurs only for J > Jps > 3b, 2- The effective binding 
energy between pairs can be quantitatively defined in the low doping limit by 



Exact Diagonalisations of 2D Correlated Fermions 



Fig. 4. Binding energy A4 vs J for clusters of 16, 20 and 26 sites (from Ref. 18). 
The 26-sites cluster hole-hole binding energy A2 is also indicated by open stars. 




J/t 

A4 = E^h + Eoh — 2£ , 2h where E„h is the ground-state (GS) total energy for a 
system with Nh = n. GS energies of 4-holes can be calculated for the three possible 

Table 2 

Binding energies A2 and A4 calculated on a 26-site cluster. 



J/t 


0.1 


0.3 


0.5 


1 


A 2 


0.05953 


-0.15496 


-0.34020 


-0.73861 


A 4 


0.38198 


0.10378 


-0.03776 


-0.22412 



orbital symmetries of the wave function (with zero total momentum) and the cor- 
responding data for A4 are shown in Fig. 3. The preformed pairs bind (ie A4 < 0) 
when J exceeds some critical value which depends on the orbital symmetry. The 
physical character of each wave function can be identified from the magnetic en- 
ergy lost in the large J limit (see Fig. 3). A first order transition is seen around 
J ~ 2.7 between a state with lines of holes to a state with holes forming a real 
droplet. The first order character of the transition in fact explains the rapid change 
in the hole-hole correlations observed recently by approximate methods. 17 In order 
to estimate the critical J A4 has been evaluated on various lattices of size up to 
N=26 and the data are displayed in Fig. 4 and Table 2. Note that the critical value 
Jb,4 at which A4 changes sign depends weakly on the system size while the slope 
|9A4/9J| at this point decreases for increasing size. Fig. 4 then strongly suggests 
that if J/t £ [J/t\ B ,2, J Mb, 4] with J/t\ B<2 ~ 0.16 and J/t\ BA ~ 0.45 individual 
pairs exist without forming larger clusters of holes. Note however that J/t\ By 2 for 
pair formation could in fact be smaller as suggested by the observation that the 
dynamical response of the d-wave pair operator exhibits on small clusters a sharp 
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(5-peak of weight Z^n cx J/t down to very small values of J/t 12 as shown in Fig. 
1. The formation of droplets should not be confused with real phase separation. 
More insight could be obtained from the calculation of the inverse compressibility 
k _1 defined by k _1 = d 2 (Ecjs j 'N) j 'd(rih) 2 . We expect k _1 to drop discontinuously 
to zero in the phase separated regime. A non-uniform phase with domains walls 
could also compete with phase separation. A careful finite size scaling is necessary 
to investigate such possibilities and the issue of phase separation remains certainly 
an open question. 
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